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Highly compressed nuclear matter created in relativistic heavy collisions 

is to large extent governed by local non-equilibrium. As an idealized scenario 
colliding nuclear matter configurations are studied within both, relativistic 
mean field theory and using more realistic in-medium interactions based on 
the Dirac-Brueckner T-matrix. The equation of state in anisotropic mat¬ 
ter is thereby governed by two competing effects: The enlarged phase space 
volume in colliding matter tends to soften the internal potential energy of 
the subsystems whereas the relative motion of the two currents leads to a 
strong additional repulsion in the system. An effective EOS constructed for 
anisotropic momentum configurations shows a significant net softening com¬ 
pared to ground state nuclear matter. This effect is found to be to large 
extend independent on the particular choice of the nuclear interaction. A 
critical discussion of standard transport approaches with respect to the con¬ 
sidered non-equilibrium effects is given. 


I. INTRODUCTION 

One major goal of relativistic heavy ion physics is to explore the behavior of the nuclear 
equation-of-state (EOS) far away from saturation, i.e. at high densities and non-zero tem¬ 
perature. Over the last three decades a large variety of observables has been investigated 
both, from the experimental and theoretical side motivated by the search for the nuclear 
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EOS. The collective particle flow is thereby intimately connected to the dynamics during the 
compressed high density phase of such reactions E.g., the elliptic flow which develops 
in the early compression phase is thought to be a suitable observable to extract information 
on the EOS p. But also the production of strange particles is a good probe to study dense 
matter |^. Recent precision measurements of the production at SIS energies strongly 
support the scenario of a soft EOS 

However, the difficulty to extract information on the EOS e(p) from heavy ion reactions 
lies in the fact that the system is out of equilibrium. At intermediate energies the relaxation 
time needed to equilibrate coincides more or less with the high density phase of the reaction. 
Hence, non-equilibrium effects are present over the compression phase where one essentially 
intends to study the EOS at supra-normal densities. Experimental evidence for incomplete 
equilibration even in central collisions has recently been reported in ||^. In the early stages 
of the reaction, the phase space distribution in the participant zone is similar to two fluids 
of counter-streaming matter. The local momentum distributions of such colliding nuclear 
matter, called CNM in the following, are given by two Fermi ellipsoids, i.e. two boosted 
Fermi spheres separated by a relative velocity |^. In the course of the reaction the mid¬ 
rapidity region is more and more populated due to binary collisions and the system is heated 
up. The originally cold and sharp momentum ellipsoids become diffuse and merge together. 
The energy density depends in the course of the reaction therefore on both, the local relative 
velocity and the local temperature: 

e(p, Urei, T = 0) —> e(p, v^ei, T) —> e(p, n^ei = 0, T) . 

In order to learn something about the ground state EOS one has to seperate the compression 
energy from the contributions of the relative motion and the thermal energy 

"^rel; T) *^comp T ^rel T ^^therm 

This is generally a difficult task since the various contributions are not independent from 
each other. Moreover, the conversion of kinetic energy into compression and thermal energy 
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goes along with the production of entropy. While the total energy is conserved, the free 
energy tends to minimize. This temporal evolution of a heavy ion reaction is described by 
microscopic transport models like BUU or QMD |^. Such simulations ||TT| , p!^ , p!3| show 
also that the local phase space in the overlapping overlapping zone of two interpenetrating 
nuclei can to large extent be idealized by colliding nuclear matter. CNM provides a smooth 
transition to a hnal, at least in the central cell, fully equilibrated spherical momentum 
conhguration. 

The aim of the present work is to investigate in more detail how the compression energy, 
related to the ground state EOS, can be extracted in kinetic non-equilibrium situations and 
which consequences arise if one tries to do this by the use of dynamical transport simulations. 
In the limit of vanishing temperature the compression energy is just the free energy of the 
system ecomp = ^{Q,Vrei,T = 0) — £^rei- To obtain an effective EOS in CNM one has to 
subtract the energy of the relative motion. Such a treatment neglects temperature effects 
but it allows to study non-equilibrium features for simplihed conhgurations. These reflect in 
hrst order the phase space evolution in the initial and the high density phase of the reaction 
where termperatures are still moderate. Therefore the inhnite counter-streaming two-fluid 
scenario has been a subject of a variety of theoretical investigations 1,00,00,130,0. 

In transport calculations the nuclear mean held is usually based on phenomenological 
parameterizations P,00,0. Such parameterizations allow different extrapolations to high 
densities, summarized by referring to a ’hard’ or a ’soft’ equation-of-state |l[], which can 
be tested in heavy ion collisions. A more microscopic approach is to start from free NN 
interactions, determine the correlated two-body interaction in nuclear matter and use this 
interaction in heavy ion reactions. Doing so, in each time step in the course of the reaction 
the mean held, respectively the EOS, would have to be derived from an integration of the two- 
body interaction over the actual momentum distribution. This determines then both, the 
density and the momentum dependence of the mean held P, |23| , 0 , |25[] . In phenomenological 
approaches, even when momentum dependent interactions are included, part of the arising 
non-equilibrium features is missing. Also with phenomenological two-body interactions the 
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complete single particle potential should consistently be derived for the actual momentum 
conhguration. In standard transport calculation this is, however, not done: the single 
particle potential contains a purely density dependent part which, as a density functional, 
has been determined for equilibrium matter. The momentum dependent part, on the other 
hand, accounts in lowest order for anisotropy and temperature effects of the phase space 
distribution. However, the density dependent part is that one which is most directly related 
to the ground state EOS. 

To obtain a model independent picture the present investigations are based on different 
choices of the nuclear forces. We use the framework of relativistic mean held theory, i.e. 
the Walecka model p6| , ^ and the microscopic relativistic Dirac-Brueckner (DB) approach 
p8| , |^ . The mean held approach has the advantage that most integrals can be evaluated 
analytically which makes the discussion of typical features of CNM more transparent. The 
DB model, on the other side, is more realistic. It turned out to be quite successful in the 
description of nuclear matter bulk properties. DB forces have been extensively tested at SIS 
energies below 1 A.GeV, and a general agreement with corresponding how data has been ob¬ 
served P^ , P5| , |5D| . Other studies of CNM have been performed in a non-relativistic framework 
using either Skyrme interactions |TB|,p!3| or Brueckner G-matrix forces ITIIT^. Relativisti- 


cally colliding matter has been mainly treated within the Walecka model ||8|Jlq J T7| J l^ . In P] 
a hrst attempt was made to apply relativistic DB forces. In some cases also hnite temper¬ 
ature ehects have been studied pO| , p!^ . In T dependent G-matrix forces were used in 
transport simulations but the inhuence of the explicit T dependence of the mean held on 
the reaction dynamics was found to be moderate. 

The paper is organized as follows: In Sec.2 we discuss the general structure of the 
relativistic energy-momentum tensor in colliding nuclear matter within both, the DB and the 
mean held (ME) approach. In Sec.3 we derive an ehective EOS for the CNM conhguration. 
In Sec.4 implications for heavy ion reactions are discussed, in particular with respect to 
transport calculations. In Sec.5 we give some short remarks concerning the treatment in 
non-relativistic approaches and hnally conclude in Sec.6. 
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II. ENERGY-MOMENTUM TENSOR IN COLLIDING NUCLEAR MATTER 


The phase space distribution of colliding nuclear matter (CNM) 

©12 = ©1 + ©2 ~ ©1 ■ ©2 ( 1 ) 

is composed by the momentum distributions of two counter-streaming matters ©* = 
©(/i*(fci?J — klu'{), i.e. two boosted Fermi ellipsoids. © is the step function, kpi are the 
Fermi momenta and u\ = ( 7 j, 7 iUj) are the streaming velocities of the two subsystem cur¬ 
rents. The last term in eq. (|l|) ensures that the Pauli principle is fulhlled for small velocities 
where the two ellipsoids might overlap. The total baryon density has thereby to be restored. 
Details of this procedure which can be performed by a covariant construction can be found 
in §. 


A. DB approximation 

The interaction of nucleons with the surrounding medium is expressed in terms of a 
self energy which dresses the particles, respectively the in-medium propagator. Analogous 
to spin-isospin saturated nuclear matter at rest, the self energy in colliding matter can be 
decomposed into scalar and vector components expressed covariantly by 



(2) 

= So(fc)M^i2 + T.v{k)^i2ku ■ 

(3) 


Here, the Sj with i = S', 0, n are scalar functions. 

The four-velocity Ui 2 is the streaming velocity of the combined system derived from the 
total current ji 2 = P 12 M 12 . Since Q 12 is the invariant baryon density in the c.m. frame of 
the two currents, gu = \fjl 2 =< 1 >12 \c.m., it is natural to work in that frame, i.e. the 
frame where the sum of the spatial components of the two counter-streaming currents vanish 
ji2 = 0 . With the help of the projector perpendicular to U12 

AfJ = g"" - (4) 
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eqs. (1^) are manifestly covariant. Inversely, the self-energy components are obtained from 
covariant projections 


Es = jtr[El 

E„ = ^tr |tif27„E] 




tr 


(5) 

( 6 ) 
(7) 


One should note that in this form the projections are identical to nuclear matter at rest 
28| except for the fact that all quantities are determined in the colliding system. Hence the 


effective mass and the kinetic momenta are given by 

M* = M + T.s{k) 

K = K + 

E* = 


( 8 ) 

(9) 

( 10 ) 


The self-energy has now to be derived from some model, e.g. from an effective two-body 
interaction like the in-medium T-matrix. In this general case the self-energy components 
(3-|l^) are momentum dependent. We also consider the mean held approach of the Walecka 
model [^] where the self-energies are not explicitely momentum dependent. This has the 
advantage that all integrals can be solved analytically. In the DB approach it is useful to 
introduce an additional conhguration averaged and momentum independent effective mass. 
In the c.m. frame (ji 2 = 0) one can thus get rid of the contribution by rescaling the 
helds in the same way as it is done in nuclear matter p8| , p9| . The spatial component of the 
self-energy (|^) has in the c.m. frame the form S = E„(fc)k. In NM has been found to be 
small and nearly constant as a function of momentum. In mean held theory this contribution 
vanishes identically. Here it can be absorbed into a rescaled ehective mass 

M -\- Eg 


M* = 


( 11 ) 


1 + 

and the rescaled momentum k* = k*/(l -|- E^) = k|c,m.- The conhguration-averaged scalar 
and vector self-energy components are dehned as 


6 








Es = < Ss{k)M’/E' >u Iqs„ 

E„ = < E„(i)/£* >12 / < l/B* >12 . 


( 12 ) 

(13) 


Here 


< X >12 — 


K 


(27r) 


<f\iX{k)e^2{k-x) 


(14) 


denotes the summation over all occupied states and E* is defined as in (p!0D, however, with 
the rescaled quantities M* and k*. In spin-isospin saturated nuclear matter the phase space 


occupancy factor in (|T^ equals k = 4. The set of parameters which determines the colliding 
configuration is denoted by y = kF 2 ,ui,U 2 , M*}. Thus the effective mass M*, the 
scalar density gsi 2 =< M*/E* >12 and the configuration (|lD itself are coupled by non-linear 
equations. The Dirac part of the in-medium propagator has the same form as in nuclear 
matter ra,pl|l 


Guik) = + M*]27n6{k*^ - M*^)2e{k*)eu{k; x) ■ 


(16) 


It should be noticed that the usage of an averaged effective mass in the propagator is a 
standard procedure of the DB approach |2^,^ (reference spectrum approximation). It is 
based on the observation that M* is only weakly varying with momentum inside the Fermi 
sea and it simplifies the solution of the Bethe-Salpeter equation considerably. To account 
for the full momentum dependence of the effective mass in the two-body propagator of the 
BS-equation is presently beyond the scope of DB calculations. In CNM this approximation 
is more severe due to large relative momenta which can occur between the two currents. 
However, since the effective mass shows in DB calculations only a rather moderate momen¬ 
tum dependence at large momenta above the Fermi surface this approximation is still 
justified. The self-energy components follow now from covariant projections of the T-matrix 
elements. For on-shell scattering of positive energy states the T-matrix can generally be 
decomposed into five Lorentz invariant amplitudes (scalar, vector, tensor, axial vector and 
pseudo-scalar or pseudo-vector, respectively) where direct (Hartree) and exchange 
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(Fock) amplitudes are related by a Fierz transformation If the amplitudes are already 
anti-symmetrized the ’direct’ representation for the T-matrix is sufficient 


S(/c) = -iJ ^^tr[G'i 2 (g) T{k,q-,x)] 


K 




M*Ts{k,q;x)+ A*Ty{k,q;x) © 12 ( 9 ; x) • 


( 16 ) 


(27r)3 J E*{q) ^ 

The scalar Ts and the vector Ty amplitude are Lorentz invariant functions which depend 
(besides the parametric dependence on the conhguration expressed by x) the c.m. mo¬ 
mentum, the c.m. relative momentum and the c.m. scattering angle between k and q in 
a Lorentz invariant manner which can be converted to a dependence on the Mandelstam 
variables 

The energy momentum tensor in CNM is given by 

d^k , , f r d^k dS 


-r'‘ve, 2 (A^;x) 


(2i)< 

K r d^k 


(27r)4 (27r)4 


tr [Gu{k) tr [Gi 2 {q)T{k, g; y)]] 


(27r)3 7 E*(k) 

K 




2 (27r)6 J E*{k)E*{q) ^ 

With eqs. (|8-W) we express in terms of the helds 


M*^Ts{Kq-,x) + k*J*^Ty{k,q-,x) 0i2(fc; x)0i2(g; x) • (17) 


=< k-"‘k-'‘'/E- >12 -I/"" - Ig'”’ (Eg 


( 18 ) 


The scalar contribution (^) and the terms containing the vector held are given by 


yfkE k*^Y.%k)lE* >12 . ( 19 ) 

The total scalar density and the total baryon current are given by 

r d^k 

es„ = -i] p^tr [G,2(4)] =< M7£- >,2 (20) 

r d^k 

Jm12 = -^ j (^^’^ b^^Gl 2 {k)] =< kyE* >12 (21) 
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B. Mean field approximation 


In order to examine the strnctnre of the energy-momentnm tensor ( 0 ) it is instructive 
to consider first the mean field approximation. To do so we use the Walecka model (QHD-I) 
26| . In QHD-I the T-matrix amplitudes are replaced by corresponding coupling constants 
for a scalar a and a vector co meson field, T 5 1 —^ T^ = gl/rn\ and Ty '—»• Ty = 
and thus the double integrals in (P!7D decouple. One could also use more sophisticated model 
such as density dependent hadron field theory [^] where r 5 ^y(p) are density dependent 
coupling functions. Such a density dependence can e.g. be taken from effective field theories 
|3^ or again from the DB T-matrix [^. In the latter case this provides the mean field 
approximation to the DB problem outlined above. In the following we consider the simpler 
case of constant couplings r 5 y but the considerations can easily be generalized. Since 
the self-energies are now momentum independent one has k* = k* and M* = M*. The 
energy-momentum tensor reads 


=< k*^k*’^/E* >12 -TvJ^2Ji2 - - ryJ^2JXl2} ■ 


( 22 ) 


In nuclear matter one obtains the kinetic energy and pressure densities as 

Ckin = < ^*(k) >= ^E^gikp) + ^M*gs (23) 

Pkin = - < k*‘^/E* >= -Epg^kp) — -M*gs . (24) 

The corresponding expressions in CNM are similar and transparent as long as the momentum 
distributions of the two currents do not overlap. Otherwise the contributions arising from 
the Pauli correction in (|I]) complicate the expressions considerably. Thus we first discuss 
the case where the relative velocity of the two currents is large enough to seperate the 
ellipsoids, i.e. 0i(A:, mi)02(A:,M 2 ) = 0. The integrals for the kinetic energy density can be 
solved analytically and yield [§] 

< k-vyE- >„= y" I Urns., + E jEsMkF.) I - E (25) 

1^ i=l,2^ I *= 1.2 
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Here Qoikpi) and EF^ = ^kp. + are the rest densities and the corresponding Fermi 
energies of the snbsystems and is the projector (||) on the snbsystem velocities. To 
fnrther simplify (|^) we consider symmetric conhgnrations = kF,ui ^2 = ±m). The 

total (c.m.) cnrrent is then given by j^i2 = + jfi 2 = (27(m)^o(^f), 0 ). Analogonsly the 

total scalar density reads 


Qsi 2 — QSi + QS2 — 


2k 


M* 


kpEp - M*2ln 


(kpE Ep \ 
M* j 


(26) 


with M* = M — TsQsi 2 - Energy density e = transverse pressnre p± = and 

longitndinal pressnre p\\ = are given by 

= 2 [^Epgoikp) + qs{M*) + Epp^ikp) ( 7 " - 1 ) 

+ ^ {Tv {2jgo{kp)f + Ts {2gs{M*)f} ( 27 ) 

p_L = 2 gs{M*) 

+ ^ {Tv {2^go{kF)f - Ts {2gs{M*)f} ( 28 ) 

P\\ = 2 |-Ei?^?o(^F) — Qs{M*) + Epgo^kp) (7^ — 1) 

+ ^ {Tv {2^go{kF)f - Ts {2gs{M*)f} ( 29 ) 


To compare CNM with NM it is instrnctive to rewrite Eqs. (00)- The separation of 
projectile and target nncleons in momentnm space increases the phase space volnme in CNM. 
The boosted ellipsoids are elongated in longitndinal direction bnt in transverse direction the 
Fermi momenta kp of the snbsystems (which are hxed by the cnrrent rest densities goikp)) 
are not enhanced. This featnre can be acconnted for if we integrate over one Fermi sphere 
at rest 0(fc; kp^u = 0), however, with a donbled phase space occnpancy factor 2k 


< X >2k — 


2k 

(27r)3 


(f'kX(k)Q{k; kp, u = 0) 


(30) 


This description yields directly the CNM scalar density (p6|) 


M* 

0Si2 =< >2 k 


( 31 ) 
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and 2gQ{kF) =< 1 > 2 k- One can now rewrite Eqs. (^^-|2^) separating thereby the static 
parts from those contribntions which depend on the streaming velocity. Making nse of (|30D 
the static parts of the energy and pressnre densities are expressed as in gronnd state matter, 
however, with k = 8: 


T“> = T“>|2, + + Tv i,(.kF) 

Pi = P±|2k + TV slikp) 

P|| = P||| 2 « + 27^u^ lEppoikp) + Tv Po(4f) 


(32) 

(33) 

(34) 


Eqs. (p^ M) are exact as long as the ellipsoids do not overlap. It can be directly seen 
from there that CNM is characterized by two effects which act in opposite direction: The 
hrst one is the separation of the nncleons belonging to the different cnrrents in phase space. 
This effect acts like an additional degree of freedom. Evidently it rednces the pressnre 
of the Fermi motion. Since the scalar part of the potential energy is not affected by the 
spread in momentnm space, besides the fact that the contribntions of projectile and target 
are snperposed, it deepens the scalar potential. However, the vector repnlsion gives rise 
to an additional velocity dependent pressnre in both, transverse and longitndinal direction 
(terms proportional to Ty). The elongation of the boosted Fermi ellipsoids in longitndinal 
direction enhances the baryon vector density which leads to an additional kinetic pressnre in 
longitndinal direction ( 27 ^u^£'i 7 ’Po(^F))- Both these contribntions are related to the kinetic 
energy of the relative motion in the system and enhance the energy density. 

The comparison to eqnilibrated NM shonld be performed at identical rest densities, i.e. 
the two cnrrents with their rest densities Po(^f) shonld be compared to NM at the same 
total density 


QNMikptot) — 2 po (^ f ) ; kpt^t — 23kF 


(35) 


By this procednre one prevents a distortion of the comparison by the pnrely kinematical 
enhancement of the boosted snbsystem densities in the c.m. frame which does not affect 
the transverse degrees of freedom. A comparison of CNM and NM at identical total baryon 
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densities would mean to interpret the Lorentz contraction in the anisotropic system as a com¬ 


pression when related to the isotropic conhguration. Condition (p5D is also necessary for a 
meaningful comparison of CNM at different streaming velocities. Doing so, one can estimate 
the phase space effects comparing the contribution T^^\ 2 K,{kF) in (|3^ ) with T^'^\FiM{kFtot) 
NM. Schematically, the separation of phase space in the colliding system is illustrated in 

Fig. I. 



FIG. 1. Schematic representation of the phase space in nuclear matter (a), colliding nuclear mat¬ 


ter (b), and in colliding nuclear matter as experienced by transverse degrees of freedom, respectively 


by the scalar, density dependent part of the interaction (c). 


It is clear that the introduction of a new degree of freedom, e.g. going from neutron 
matter to isospin saturated matter, strongly increases the binding energy. In contrast to the 
scalar held, where the colliding conhguration can be exactly treated by the introduction of 
an additional degree of freedom (Fig.|l|ai— ^-c), the vector held is enhanced by the spread in 
momentum space (Fig.|l|b) and the elongation of the boosted ellipsoids (not shown in Fig.|Ilb). 
However, the present picture is more general. Besides the fact that in CNM the internal 
Fermi pressure of the subsystems is reduced also the contributions of the internal potential 
energy of the subsystems are modihed due to the smaller Fermi momenta. The spread in 
momentum space, on the other hand, determines the interaction between the two currents 
and leads to additional repulsion in the combined system. These are the two essential 
features of colliding matter and - as will be seen later on - their implications on the EOS 
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are relatively independent on the choice of the nuclear forces. 

The above considerations (Eqs. (p6| - p^) ) required that the ellipsoids were completely 
separated. The overlapping case which contains the limit Urei —^ 0, i.e. equilibrated NM 
is more complicated. This involves the Pauli correction and the construction of the CNM 
conhguration without double occupancy [§. The principle of this procedure is to keep the 
streaming velocities Ui hxed and to enlarge the corresponding Fermi momenta {kp —^ kp) 
thus that the condition 


joi2{kp) — 


K 


(fk 


(27r)3 . 

= ‘^l{u)Qo{.kp) 


0i(/c; kp, u) + 02 (A;; kp, —u) — 0i(fc; kp, M) 02 (fc; kp, —u) 


(36) 


is fulhlled. Eq. (p6f) is valid for the symmetric case but can easily be extended to asymmet¬ 
ric conhgurations 0]. In the mean held approximation the integrals over the intersection 


ellipsoids can be worked out analytically which has been done by L. Sehn [^]. For the 
energy density one ends up with 


T^^ = -M\kp) i2gs{kp)-6gs) 

+ X \^v l‘^{‘2go{kp) — SgoY + Ts {2gs(kp) — (5p5)^| 


+ 


K 

Stt^ 


1 --M*Hn 

jl + u Ep — kp 

1 2 

1 

1 

s 

* 

1_ 


kpEpM*"^ + —'y‘^Ep{Ep — ukp)^ — — ( - 

3m 3m \ 7 


-j^M*^(uEp — kp){Ep + ukp) 


(37) 


The Pauli correction terms 6gs and 6go (originating from the 0i02 overlap integrals) can 
be found in [Q. In the limit m —> 0 the CNM conhguration (|I]) approaches smoothly 
the corresponding NM conhguration Q{kp^^^) with kp —^ kp^^^, ^Qsikp) —>■ gs{kp^^^) and 
Sg^ikp) —>• go{kp^^^). The same holds for the energy momentum tensor, i.e. T^^''{kp^u —>• 
0) —>■ T^'^{kp^^^,u = 0). The enhancement factor which accounts in (j^-^) ehectively for 
the enlarged phase space volume approaches unity, i.e. 2k —> k. 
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III. THE EFFECTIVE EOS 


In colliding nuclear matter the energy per particle is defined in the usual way 


^12(^12,m) = T°7pi2-M (38) 

with Qi 2 = \fj \2 the invariant baryon density in the c.m. frame of the two currents. 0 
However, a meaningful discussion of non-equilibrium effects with respect to the ground 
state EOS should be based on the binding energy per particle. In the temperature zero case 
this corresponds to the free energy of the system where the contribution from the relative 
motion of the two currents has been subtracted. The corresponding effective EOS in colliding 
matter is then also directly linked to the hydrodynamical picture |]I|. The kinetic energy of 
a nucleon inside the medium is given by = E*{k) — M*{k). Analogously, we identify 
the kinetic energy density of the relative motion in (|T^ with 

4in=<fc*°(fc*°-M*)/E*>i2=<E*(fc)-M*(7 >12 . (39) 


The effective mass M* accounts in (p^ for the fact that the two currents are interacting 
by momentum dependent forces. Such a momentum dependence occurs in the relativistic 
treatment already on the mean field level where the self-energies do not explicitely depend 
on momentum. The magnitude of M* characterizes thereby the strength of the momentum 
dependence of the optical potential. On the other hand side, different parameterizations 
of QHD including additional non-linear self-interaction terms for the scalar a field la can 
yield identical EOSs but differ in their momentum dependence. This effect is taken into 


account in (|39|) . Moreover, the kinetic energy equals E* — M* ~ fhe non-relativistic 


^According (|^) the c.m. density is related to the subsystem rest densities by ^>12 = 2j{u)Qo{kF). 
To avoid the stretching of the density scale by the Lorentz contraction we will in the following 
represent the energy density as a function of the subsystem rest densities 2QQ{kF). Physically this 
means to remove the Lorentz contraction in the representation of the free energy. 
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limit. In the non-relativistic limit M* plays the role of the Landau mass which is used to 
parameterize the optical potential. 

The energy of the relative motion of two interacting currents is more subtle to dehne. 
The present description is inspired by the fact that in the non-relativistic case this should be 
the energy which is necessary to shift the two Fermi spheres in their c.m. frame against each 
other. In the relativistic framework this corresponds to the difference of the kinetic energy 
in the combined system and the corresponding NM conhguration at twice the subsystem 
rest density 

S,ex{kF,u)=<E*-M*>^2-<E*-M*>^=^ . (40) 

By dehnition f’j.ei contains the kinetic energy which arises from the separation of the two 
distributions, respecting thereby the Pauli principle, and accounts for the interaction be¬ 
tween the two currents by the presence of the effective mass. The vector density is not 
conserved when boosting the two distributions against each other. The Lorentz contraction 
contributes therefore to the energy of the relative motion. Similar as the c.m. energy density 
(e = r°°|c.m. = Eq. (^D]) can be formulated in an invariant manner 

Srelikp^u) = Uu^ < - M*6;^)/E* >12 Uuu 

• ( 41 ) 

Now one obtains the binding energy per particle 

rpOO _ C 

E'f^^ikF, u) = -- M (42) 

Q12 

as a function of the total c.m. density pi 2 and the c.m. streaming velocities. In the following 
we will only consider symmetric conhgurations j = Eu). The basic features 

of the colliding conhguration are easier discussed at the mean held level where formulas have 
been evaluated analytically and thus we consider this case hrst. 
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A. Mean field approximation 


Here we mainly use the QHD-I version of the Walecka model with only two free param¬ 
eters adjusted to the nuclear matter bulk properties. The model yields a binding energy of 
-15.75 MeV, a relatively large saturation density = 1-42 fm~^), a small effective mass 
of M*/M=0.556 and a large compression modulus of K=540 MeV. 



u/c 

FIG. 2. Effective mass in colliding nuclear matter (CNM) at subsystem densities go = Qsut as a 
function of the streaming velocity (solid) compared to the corresponding effective mass in nuclear 
matter obtained at with a doubled phase space factor (2k, dashed). Results are obtained with 
the Walecka model QHD-I. 


To illustrate the influence of the two-current geometry on the scalar part of the inter¬ 
action in Fig. H the effective mass in colliding nuclear matter at subsystem densities Psat is 
shown as a function of the streaming velocity. At zero velocity the nuclear matter result at 
twice saturation density (p = 2psat) is obtained. With increasing streaming velocities the 
effective mass drops and reaches around u/c = 0.7 (which corresponds to i7iab=l-8 GeV) the 
asymptotic value of the non-overlapping configuration given by expression (^^. A further 
increase of u has no influence on the effective mass since the scalar densities can be evaluated 
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in the local rest frames of the currents. At the mean held level M* is exclusively determined 
by the shape of the momentum distribution, but once the ellipsoids are separated, the ef¬ 
fective mass stays constant. The asymptotics are given by the corresponding effective mass 
in nuclear matter at go with the enlarged phase space factor {2k,). Hence the two-ellipsoid 
geometry enhances the attractive scalar held by about 30 MeV. However, in the total energy 
per particle this enhancement is completely compensated by the additional repulsive 
contributions in the counter-streaming system. As discussed above, Eqs. (|3^|34|) , an addi¬ 
tional longitudinal Fermi pressure occurs due to the boosted vector densities. Secondly, the 
vector held gives rise to a strong repulsive force between the two currents. 


200 
150 

> 100 

Qi 

^ 50 
0 

-50 

0 12 3 12 3 

P/Psat 

FIG. 3. Energy per particle E 12 in colliding nuclear matter as a function of the subsystem 
density 2po/^>sat at different streaming velocities tt/c=0/0.2/0.4/0.6. Results are obtained with the 
Walecka model QHD-I (left) and the non-linear Walecka model NL2 (right). In both cases densities 
are normalized to the saturation density ^sat = 0.193 fm“^ of QHD-I. 

As can be seen from Fig. the additional scalar attraction in the total energy per particle 
(^) is then completely hidden by the much stronger repulsive contributions. Already at 
very low streaming velocities the system appears to be unbound. One should, however, be 
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aware that QHD-I is a model with an unrealistically strong repulsive character at higher 
momenta, expressed by its small effective mass and the large vector held. As a result QHD-I 
strongly over predicts the empirical optical nucleon-nucleus potential and the transverse 
nucleon how in heavy ion collisions P,p6| . To estimate the model dependence inherent in the 
present considerations in Fig. ^we compare also to a softer version, the non-linear Walecka 
model NL2 |^, with bulk properties oi E = —15.75 MeV, = 1.29 fm“\ K=237 MeV. 
A large ehective mass M*/M=0.80, respectively a weaker momentum dependence leads to 
a more realistic description of the dynamics in heavy ion reactions [p|,^. Although less 


repulsive at large momenta the NL2 parameterization yields qualitatively the same velocity 
dependence of the total energy per particle as QHD-I. 



FIG. 4. Transverse pressure p±_ in colliding nuclear matter as a function of the volume per 
particle at different streaming velocities. Results are obtained with the Walecka model QHD-I 
(left) and the non-linear Walecka model NL2 (right). 

In Fig. ^ the transverse pressure pp = is shown as a function of the volume 

per particle V/A = 1/Qi 2 - Interestingly, the transverse pressure is reduced at low streaming 
velocities. Although the momentum distribution of the two currents is still quite compact at 
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u/c = 0.2 the reduction of the Fermi pressure and the reduction of the effective mass by the 
enlarged phase space volume are sufficient to produce this softening in transverse direction. 
However, with increasing velocity the vector repulsion starts to dominate. Qualitatively this 
behavior is model independent, i.e. it occurs in QHD-I and the softer NL2 model. Only at 
the highest velocity shown in Fig. ^ {u/c = 0.8) the ellipsoids are well separated, i.e. eq. 
(^) is exactly valid. At smaller u the Pauli corrections are taken into account. 

In order to investigate the subtle features of the colliding system we consider now the 
effective binding energy per particle With eq. (^) the energy of the relative motion 

can easily be evaluated in mean held approximation (no overlap): 


EFh^-T)- - lM‘ + 


Sj-el 2pg(/Cp) 

+ ^ [M*gs 12 - M*{kF^^Jgs{kFtot)] 


(43) 


The meaning of (^3l) becomes more transparent in the low density limit. Using 


0s — 00 
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(44) 


one obtains after an expansion in u/c 
£^rel ~ 2 Po ( ) 


M* , 3 e 
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2 5 2M* 2M* 


^Ftot 


+ 0{u^ 


(45) 


5 2M*{kF,J_ 

Now it becomes evident that eq. (]45|) contains the kinetic energy of two freely streaming 
gases — 7 M*po- The hrst term is just the mean kinetic energy of a comoving nucleon, the 
second term is the kinetic energy of the internal Fermi motion in the local rest frame and the 
third contribution is due to the elongation of the Fermi ellipsoids in longitudinal direction 
by the boosts. Subtracted is the kinetic energy of the Fermi motion of the rest system at 
equal rest densities. Hence the dehnition (^31) respects the Pauli principle. An alternative 
choice would be to subtract simply the kinetic energy of the superimposed two spheres as 
illustrated by Fig. |^c which would, however, violate the Pauli principle. To subtract the 
kinetic energy of the rest system at identical, Lorentz enhanced vector densities, on the 
other hand side, would introduce an additional unreasonable velocity dependence into the 
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last term of eq. (^51) . Thus eq. (^) is the natural dehnition of the relativistic kinetic energy 
of the relative motion and contains the correct non-relativistic limit. 



of the subsystem density 2po/£'sat at different streaming velocities. Results are obtained 


with the 


Walecka model QHD-I (left) and the non-linear Walecka model NL2 (right). In both cases 


densities 


are normalized to the saturation density ^sat of QHD-I. 


Now one is able to construct an effective EOS which relates colliding to ground state mat¬ 
ter. In Fig. 1^ the corresponding effective equations of state are shown. First of all it is seen 
that the two-Fermi-ellipsoid geometry leads to a softening of the effective EOS at moderate 
streaming velocities. This can be understood by the reduced Fermi pressure in transverse 
direction and the enlarged scalar attraction as compared to NM. The vector potential is 
not affected by the subtraction scheme for £^rei and contributes fully to Thus, the 

enhancement of the vector repulsion acts in opposite direction and makes the EOS harder 
at large densities and/or high streaming velocities. Here the model dependences become 
more significant. At u/c = 0.8 the effective QHD-I EOS lies well above the corresponding 
ground state result whereas for NL2 it is still softer at low densities. This reflects nicely the 
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weaker momentum dependence and the weaker repulsion of NL2. In summary the softening 
of the “’effective EOS” ’ by about 10 MeV due to the anisotropic phase space distributions 
in CNM is hardly seen on the scale of the total energy per particle but on the scale of an 
effective binding energy. Such small effects are, however, the features one is after in heavy 
ion experiments at intermediate energies. 


B. DB approximation 


The mean held used in the present section is based on recent DB calculations with 

nuclear matter saturation properties (Bonn A potential) of psat = 0.185 fm“^, = 

— 16.15 MeV, a compression modulus of A' = 230 MeV and an effective mass of M*/M = 
0.678 which lies in between the phenomenological models QHD-I and NL2. In ground 
state matter the momentum dependence of the DB self-energies is moderate above the 
Fermi surface. The solution of the full DB problem for colliding nuclear matter, i.e. the 
self-consistent solution of the BS-equation for two-Fermi-ellipsoid conhgurations, is still an 
unresolved problem. Thus we determine the effective interaction in CNM as described in [Q. 
The ground state T-matrix amplitudes are parameterized in form of averaged momentum 
and density dependent scalar and vector coupling functions 

fsik^kp) = [ S^M*Ts{k,q-,kF) 0(g;M / [ 0(g;M (46) 


F*(q) 


E*{q) 


Tvik^kp) = f ^ ^ Tv {k,q-, kp) Q^q^kp) / f ^ &{q]kp) 


E*{q) 


E‘(q) 


(47) 


The amplitudes Ts,v in eqs. ([T6|JT7|) should self-consistently be determined for colliding mat¬ 
ter conhgurations. However, in the present work they are approximated by the corresponding 
NM amplitudes read in the subsystems at the respectively transformed momenta k'. 
Details of this procedure can be found in |^. The treatment is not self-consistent on the 
level of the ehective two-body interaction, e.g. the inhuence of the anisotropic momentum 
conhguration on the Pauli operator for the intermediate states in the BS-equation is not 
taken into account, but it accounts for the correct geometry and the integration of the 
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T-matrix amplitudes over this geometry. Fig.^ displays the momentum dependence of the 
Ts,v amplitudes at various densities. At high momenta which exceed the range of the DB 
calculations we assume a l//c asymptotics for the amplitudes. 
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FIG. 6. Averaged scalar and vector T-matrix amplitudes Tg^v at various densities. The ampli¬ 
tudes are obtained from DB calculations for NM at rest. 


The corresponding equations-of-state are shown in Fig.|^ The behavior of both, the 
energy per particle E 12 and the effective binding energy per particle are similar to the 
MF models discussed above. As already expressed by the saturation properties the DB model 
is signihcantly less repulsive than QHD-I. Moreover, the explicit density and momentum 
dependence of the amplitudes Tg v leads to some visible deviations from the MF picture. 
At low streaming velocities {u = 0.2, 0.4) even the total energy Eu lies slightly below the 
ground state result. Consequently, the increase in binding energy seen in the effective EOS 
(lower part of Fig.|^) is more pronounced than in the MF approach, in particular at high 
densities. This behavior is due to the fact that the amplitudes Tg v decrease with density 
and momentum. Compared to the MF case this momentum dependence leads to a reduction 
of both, the scalar and the vector held. In particular the reduction of the vector held, which 
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increases in the MF models linearly with 7 and g, results in considerably less repulsion for 
highly anisotropic CNM conhgurations. 
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FIG. 7. EOS in nuclear matter (solid) and colliding nuclear matter as a function of the subsystem 
density 2po/^>sat at different streaming velocities. The DB model is used. The upper and lower parts 
show the total energy per particle £'12 and the effective EOS, i.e. the binding energy per particle 
where the kinetic energy of the relative motion in CNM has been subtracted, respectively. 
The streaming velocities are u=0.2 (dotted), 0.4 (dashed), 0.6 (long-dashed), 0.8 (dot-dashed). In 



both cases densities are normalized to the saturation density ^>sat the DB model. 


The momentum dependence of the T-matrix amplitudes is also reflected in the conhgu- 
ration dependence of the effective mass M* shown in Fig.^ as a function of the streaming 
velocity at hxed subsystem densities qq = Psat- As in the MF approximation, M* drops 
hrst with increasing anisotropy but does not fully reach the limit of a doubled phase space 
factor {2k). Instead, M* starts to rise again at high relative velocities. This reflects the 
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decrease of the scalar amplitude at large momenta, respectively the reduction of the scalar 
held when the momentum integration over the conhguration is performed. In summary, the 
DB model follows to large extent the behavior of the MF approach. The leading kinematical 
effects such as the separation of phase space and a dropping effective mass as a function of 
the anisotropy are present. Superimposed are the density and momentum dependence of 
the T-matrix amplitudes. The explicit momentum dependence of the effective interaction 
leads to a reduced strength of both, the scalar and the vector held at large relative momenta 
which results in an even softer ehective EOS in CNM compared to the MF picture. 



FIG. 8. Effective mass in colliding nuclear matter (CNM) at subsystem densities qq = ^psat as a 
function of the streaming velocity in the DB model. The result is compared to the corresponding 
effective mass in nuclear matter obtained at with an enlarged phase space factor (2k, dashed). 

IV. HEAVY ION REACTIONS 

In this section the connection of the previous considerations with transport calculations 
is discussed. The early and high density phase of relativistic heavy ion reactions in the 
SIS (Eiab < 2 AGeV) and AGS (Eiab < 10 AGeV) energy range are characterized by a 
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high anisotropy of the local momentum space in beam direction. This has also practical 
implications for heavy ion reactions and the interpretation of experimental observables. In 
the overlapping zone of two interpenetrating nuclei the initial conhguration is that of two 
sharp momentum ellipsoids. In the course of the reaction the mid-rapidity region is then 
more and more populated due to binary collisions and the system is heated up, i.e. the 
originally cold and sharp momentum ellipsoids become diffuse and merge together. At least 
in the central cell dually a fully equilibrated spherical configuration can be reached. This 
has e.g. been demonstrated in |]IT i|T^p!^ . 
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FIG. 9. Time evolution of the kinetic part of the energy-momentum-tensor T^in (left) and the 
baryon density (right) at the collision center in a central Au + Au reaction at 1 A.GeV. The tensor 
components parallel and transverse to the beam direction are shown separately. 

To obtain a more quantitative measure for the size and relevant time scales for phase 
space anisotropies in Fig.^ the time evolution of the kinetic part of the energy-momentum- 
tensor at the collision center of a central (b=0 fm) Au + Au reaction at 1 A.GeV is compared 
to the time evolution of the corresponding baryon density. The difference between the 
parallel ^kin = and the transverse = (T^i^+T^^^)/2 components is thereby a measure 
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for the anisotropy of the local momentum space in beam (z) direction. The quantities are 
obtained from relativistic BUU calculations where is determined from the phase space 
distribution /(x, k, t) represented by testparticles 


1 A-N 

=< k‘'‘k“'/E' >f= — . 


(48) 


The index i in Eq.(^) refers to the testparticle and 5 ^ is a covariant Gaussian |^. In the 
initial phase up to about 10 fm/c the local phase space is highly anisotropic, i.e. >> 
During the high density phase from ~ 5 -h 25 fm/c the system starts to equilibrate. At 
the end of the high density phase and during the expansion phase the two components 
of the energy-momentum tensor are almost equal which indicates that a large amount of 
equilibration has been reached. Fig.^ demonstrates also that the relaxation time to reach 
equilibrium conhgurations coincides more or less with the high density phase 

of the reaction. 

Now the question arises how accurately the BUU phase space can be reproduced by a 
sequence of colliding nuclear matter conhgurations. It is clear that the considered CNM 
conhgurations are an idealization of the real phase space. Firstly, one assumes symmetric 
conhgurations in projectile and target currents whereas the transport simulations contain 
density huctuations. Secondly, the system is heated up whereas temperature ehects are 
neglected in the cold conhgurations. To estimate the validity of this approximation the 
following procedure is applied: 

The temporal evolution of the two parameters Qo = Qq + Qq and Urei which characterize 
the anisotropic momentum space in terms of CNM are directly determined from the RBUU 
simulation, again in the central cell for Au+Au at 1 AGeV. The ehective compression density 
Po is given as the sum of the projectile and target invariant rest densities pg= \/ 
and the velocity Urei = |vrei| is dehned via the projectile and target streaming velocities 

(49) 


P T 

u — u 


Vrel = 


u^-u^ 


These values are compared to the corresponding po^'^rei parameters for idealized symmetric 
(^0 = (^ 0 ) CNM conhgurations. Here Urei is chosen in such a way that at given qq the 
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kinetic energy-momentum tensor components from the transport simulations are 


reproduced. The comparison is performed in Fig. |T^ where each symbol corresponds to 
a time step of At = 0.2 fm/c. Although not perfect, the CNM approximation is able 
to describe the time evolution in the qq — Urei parameter space quite well. Since the CNM 
parameters are adjusted to reproduce the kinetic energy, deviations in the Urei parameter can 
be regarded a measure for the difference between the real and the idealized conhgurations. 
These are due to temperature and hnite size effects as well as to asymmetries arising from 
density fluctuations. One has thereby, however, to keep in mind that a description in terms 
of a naive local density approximation would completely neglect the dependence on the Urei 
parameter. 

From the time evolution of the system in the po ~ '^^rei parameter space the following 
can be seen: Asymptotically the two nuclei start to touch at low density and high Urei. 
Then the density increases rapidly and reaches a maximum value of about 2.5 Psat- In the 
beginning Urei is large but then equilibration sets in. When maximal compression is reached 
the relative velocity is still relative high (urei/c — 0.6), i.e. the local momentum space is 
still highly anisotropic. In the expansion phase the density drops rapidly and the system 
becomes more and more isotropic. Interestingly, the system stays for some time in the 
region of highest compression before it expands, i.e. it takes some time before an isotropic 
pressure builds up. Also in it was found that only in the hnal expansion phase the 
pressure is isotropic. Thus the momentum space anisotropy is most distinct before maximal 
compression is reached, consistent with Fig. One effect which should also be noted is 
the slight increase of Urei in the hrst two fm/c. This acceleration is a typical effect due to 
the interaction between the two currents. In the interacting system the effective mass drops 
which leads to the slight acceleration before the stopping due to binary collisions sets in. 
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FIG. 10. Time evolution of the two parameters, total density qq and relative streaming velocity 
Urei, which characterize anisotropic momentum configurations. The squares denote the values 
obtained in the central cell from a RBUU calculation for a central Au + Au reaction at 1 AGeV. 
The triangles are the values determined by adjusting colliding nuclear matter configurations to the 
kinetic energy distributions of the RBUU calculation. Each symbol corresponds to a time step of 
0.2 fm/c. 


Thus one sees that the picture of two counter-streaming fluids is appropriate for a ma¬ 
jor part of the reaction, in particular when supra-normal densities are present. Certainly 
the assumption of cold fluids has to be improved since the system is heated up. However, 
with respect to the local density approximation where the Urei dependence is to large extent 
neglected in the determination of the mean held the two-huid picture provides an impor¬ 
tant step forward towards a more precise dependence of the mean held on the phase-space 
conhguration. In transport calculations like BUU or QMD or their relativistic coun¬ 
terparts RBUU and RQMD the mean held is usually determined in a sort of a local 
density approximation. In the simplest approximation where no explicit momentum depen¬ 
dent forces are considered the helds depend only on the total density. In a non-relativistic 
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approach this yields a local potential of the Skyrme type f4.p. = ctps + P{,PbV■ In the 
relativistic case the single particle potential has the form 

f/s.p.(ps, k) = k^ - VM2 + k2 = Vk*2 + - So - . (50) 

In mean held approximation, i.e. for Walecka type models, the scalar and vector helds 
depend on the local density, i.e. M* = M + Tis{pb), Sq = So(ps) and the lowest order 
momentum dependence of Ug.p., which is absent in non-relativistic treatments, originates 
from the difference between M and M*. It is worthwhile to extract the error which is 
made when the local density approximation (LDA) is applied to anisotropic momentum 
space conhgurations. Naturally, the LDA accounts only insufficiently for such a scenario. In 
transport simulation the particles are usually propagated in the center-of-mass frame of the 
colliding nuclei which corresponds to the c.m. frame of the two counter-streaming currents. 
There the baryon density is given by ps = P 12 = ‘^IQo which means that the LDA treats the 
Lorentz contraction of the nuclei like a compression. 

To estimate the resulting error in the dynamics we compare the correct single particle 
potential for CNM to that obtained by making use of the LDA. Therefore a typical nucleon 
comoving with one of the currents with its average velocity u is chosen. In the c.m. frame the 
spatial components of the total current ji 2 vanish and the vector field contains only the time 
like component. Hence the momentum of such a “test” nucleon is given by k = k* = M*'yuez- 
Since the vector helds is on the mean held level identical in both, the CNM calculation and 
using the LDA, i.e. Sq = —TvQi 2 = —hyps, diherences in the single particle potential arise 
only due the diherent values of the ehective mass. 
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FIG. 11. Single particle potential in colliding nuclear matter for a test nucleon which is comoving 
with one of the currents at subsystem rest densities of 0.5 and 1.0 ^>sat- The CNM calculation is 
compared to the local density approximation. Results are obtained with the Walecka model QHD-I 
(left) and the non-linear Walecka model NL2 (right). 

In Fig. 0 the single particle potential (|^) for snch a test nncleon moving with streaming 
velocity u is shown as a fnnction of the single particle energy Es,p, = kP — M, for the two 
subsystem densities qq = 0.5^sat and l.O^sat- In this representation both, E^.p. and the baryon 
density pb = depend on u. The effective mass and the resulting nucleon momentum are 
consistently determined for the colliding conhguration (including Pauli effects). We compare 
this to the LDA where the helds are determined as for ground state nuclear matter at ps. 
Due to the different value for M* one obtains a different test nucleon momentum at hxed u. 
As seen from Fig.|^ the LDA leads to signihcant deviations in the single particle potential. 
That the deviations are even more pronounced at lower densities (O.Spsat) comes from the fact 
that at high densities the potential is dominated by the vector held which coincides in the 
two approaches. Comparing QHD-I and NL2 illustrates the model dependence of the error 
made by the LDA. At higher energies the LDA predicts generally a too repulsive potential, 


30 








at low energies on the other hand too much attraction. The reason lies in the interplay 


between M* and k which both enter into f/g.p, and T^s.p.- For increasing streaming velocities 
at hxed subsystem rest densities the effective mass cannot drop below the asymptotic value 
reached for the completely separated ellipsoids. However, the effective mass for ground 
state matter, determined as a function of the increasing baryon vector density, decreases 
below this limit. This leads to an overprediction of the scalar attraction in CNM and a 
correspondingly smaller value of the attributed momentum k. The interplay between these 
two effects creates the differences in the mean field potential. 

A correct determination of the effective mass, respectively the scalar density is an im¬ 
portant issue, not only with respect to the nucleon potential. Also mesons like pions, kaons. 


p-mesons ect. are assumed to change their properties in a dense nuclear environment [38 


For the pseudo-scalar octet the Gell-Mann-Oakes-Renner relation connects the vacuum 
mass with the expectation value of the scalar quark condensate. Within the framework of 
QCD inspired effective models, like NJL |^, effective chiral models etc. the change 
of the scalar quark condensate in the medium is usually expressed in terms of the scalar 
baryon density. E.g. the resulting in-medium dispersion relation for kaons 


a;K±(k) = m 


-‘KN 


K 


p 

J TT 


es + vg ± v„ 


(51) 


contains an attractive scalar part 'Lkn/{IDqs and the vector potential = 3/(8/.^)j^ which 
is of different sign for kaons and antikaons. In particular kaons are primordially produced in 
the early and high density phase of the collisions which means that they are born into 
an environment which resembles to large extent colliding matter. Already at SIS energies 
but particularly at AGS energies the scalar density can be given by the asymptotic value 
(p 6 |) for the corresponding GNM configuration, characterized by high streaming velocities, 
large 7 -factors and high pB- Hence the LDA will lead to an overestimation of the scalar 
density and the corresponding mass reduction of these particles. 
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FIG. 12. Scalar density in colliding nuclear matter (CNM) at subsystem densities go = 0sat as a 
function of the laboratory energy determined in the Walecka Model QHD-I. The result is compared 
to the corresponding effective mass in nuclear matter obtained in a local density approximation to 
the colliding system. 


To estimate this effect in Fig.|^ we show the scalar density in CNM as a function of the 
laboratory energy and compare to the LDA result. The laboratory energy of the asymptoti¬ 
cally non-interacting system is given by T^iab = M(j(v) — 1) = — 1) where v is the 

laboratory velocity and u is again the c.m. streaming velocity. The subsystem rest densities 
are again chosen as go = Psat and QHD-I is used. The system is not compressed above 2po, 
i.e. the increase of pb = 2j(u)go is only due to the relative motion of the two currents. 
In the LDA the Lorentz contraction is interpreted as compression and the scalar density 
is calculated at the correspondingly enlarged Fermi momentum. Evidently this standard 
procedure yields completely different results. 

Consequences for observables, in particular collective flow observables, have been dis¬ 
cussed in [^]. The considered non-equilibrium effects reduce the transverse flow which is 
consistent with the softening of the effective EOS discussed in Sec. III. The magnitude 
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elliptic flow 


of this effect was found to be of the order as the usage of different (soft/hard) mean field 
parameterizations. 
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FIG. 13. Excitation function of the elliptic flow V 2 in Au+Au reactions at mid-rapidity. The data 


are taken from the FOPI- (triangles) and EOS/E895 (diamonds) |22] collaborations. RBUU 
calculations are performed with DB forces applied in the local density approximation (LDA) and 


determing the mean field locally in the colliding nuclear matter approximation (CNM). 


As an example we show in Fig. the excitation function of the elliptic flow V 2 for 
Au -|- Au reactions at mid-rapidity. The mean held in the corresponding RBUU calculations 


47| is based on the DB model discussed in Sec. III. In the one case the mean held is simply 


treated in the local density approximation (LDA), in the second case the potentials for CNM 
were used. Here the local parameters t) Urei(x, t) are determined from the actual 

phase space distribution /(x, k, t). Details of the calculations can be found in |2^ and [T7 . 


The two different treatments which are based on identical nuclear forces yield signihcantly 
different results for the V 2 - Since V 2 at mid-rapidity is in particular sensitive to the EOS 
at supra-normal densities these effects are here most pronounced. Using the local 

density approximation one would exclude the underlying EOS from the comparison to data 
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as too stiff. The more consistent treatment of momentnm space anisotropies leads to the 
discussed net softening of the effective EOS in the course of the reaction and restores the 
agreement with experimental data. Therefore non-equilibrium effects should be taken into 
account on the level of the effective in-medium interaction which means to determine the 
mean held used in transport calculations consistently for colliding nuclear matter. 


V. NON-RELATIVISTIC APPROACHES 


Since non-relativistic mean helds are also widely used in transport models it is worthwhile 
to shortly consider this case, in particular with respect to the structure of the single particle 
potential. For this purpose we choose a simple interaction which can serve as an illustrative 
example. The non-relativistic problem of two Galilei transformed Fermi spheres has e.g. 
been treated in []^. Using a spin-isospin averaged Skyrme force of the type 


V (k, q; k', q') =to + ^-^[{k- q)^ + (k' - q')^] + t 2 (k - q) ■ (k' - q') (52) 


the energy density in nuclear matter reads 




(63) 


Three-body forces which give rise to a quadratic density dependence of the mean held have 
been neglected in the interaction (|5^ . Thus we restrict the present discussion to the two- 
body level. The kinetic part of the energy density is given by 


with K = 4 in spin-isospin saturated matter. The local mean held and the ehective Landau 
mass are given by 


UUkp) = ^ 

OQ 


which yields the single particle potential 


1 de 
2m* dr 


( 55 ) 
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Us.p.{kF,k) = Uiocikp)+ =toQ+{ti+t2)T+ k'^{ti+t2)g . (56) 


The energy density in non-relativistic colliding matter, i.e. for two separated spheres {Mu > 
2kF) with snbsystem densities g{kF), total density QuikFtot) = ‘^Q{kF) and c.m. streaming 
velocities ±n has the form 
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(57) 


The different contribntions in Eq.(^) are easy to interpret: Like in the relativistic case the 
energy density can be decomposed into a static part e| 2 K which is determined by the donbled 
phase space volnme, and a velocity dependent part. The latter contains the kinetic energy 
of the relative motion and the potential energy which originates from the integration of the 
interaction over the momentum spread of the two spheres, i.e. the interaction between the 
two currents. With (^) one obtains the local part of the single particle potential as 




loc 12 


^£12 

dgi 2 


— Uloc 12(^f) + f7jbc 12i^F,u) . 


(58) 


/^From (^7|) and (|^) one hnds that the velocity independent part of tdoc 


Uloc 12{kF) = ^ = 2Uloc{kF) (59) 

ogi2 

equals two times the potential in nuclear matter at half the total density, i.e. it behaves 
like the scalar part of the relativistic mean held. In both cases the differences between 
the local potentials derived in NM and CNM originate from non-linear density dependence, 
(^), respectively the momentum integration over the interaction which creates this density 
dependence. In the relativistic case the M*/E* weights ([2l|j26|) play this role. Thus the 
reduced Fermi pressure affects not only the kinetic but also the internal potential energy of 
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the subsystems. Although the Skyrme force (|5^) will lead to unrealistic quantitative pre¬ 
dictions for CNM at large velocities (it was especially designed for small relative momenta) 
it serves as an instructive example for a momentum dependent two-body interaction. The 
same qualitative arguments are valid when a more realistic interaction, e.g. the Brueckner 
G-matrix, is used in (^). A stronger binding in colliding nuclear matter has e.g. been ob¬ 


served in non-relativistic G-matrix calculations for colliding matter |]I^, but without giving 
an interpretation. 

One should be aware that this type of phase space effects is not included in standard 
transport calculations for heavy ion collisions, even when momentum dependent interactions 
are used. Phenomenological potentials are usually composed by a local, density dependent 
potential and a non-local momentum dependent part . For a nucleon with momentum k 
the single particle potential derived from the phase space distribution /(q) oc ©12 (q) reads 


Gs.p, (fc_Ftot) b^loc(^Ftot) A 


(27r)3 . 

= Gloc(fcFtot) A Gnonloc(^Ftot5 


(fq r(k,q;*:f„)ei2(q) 


( 60 ) 


In transport calculations Gnonioc arises from the integration of an effective two-body inter¬ 
action over the actual momentum distribution /. Thus Gnonioc accounts properly for the 
spread in momentum space between projectile and target nucleons as well as for the reduced 
Fermi momenta inside the subsystems. The same procedure should be applied to the local 
part of the potential. However, in standard transport calculations the density dependence 
of the mean held is parameterized into Gioc as a function of the total density. This density 
dependence is decoupled from the momentum space anisotropy and differs from (pSl) . Gon- 
sequently, U\oc{kF^_^^) rehects the density dependence of equilibrated matter but should be 
calculated consistently also for anisotropic momentum space conhgurations. 


VI. SUMMARY 

The early as well as high density phase of relativistic heavy ion reactions from SIS up 
to AGS energies are to large extent governed by phase space conhgurations far from local 
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equilibrium. The phase space of the participant matter can be approximated by counter¬ 
streaming or colliding nuclear matter. In the present work we discussed implications for the 
nuclear equation of state which occurs in snch non-eqnilibrinm conhgurations. We restricted 
the discussion to the temperatnre zero case. 

The EOS in colliding matter is determined by two competing effects: The separation 
of phase space rednces the Fermi pressnre in the snbsystems compared to gronnd state 
matter. The redaction of the internal snbsystem Fermi momenta affects also the potential 
part of the energy density and leads to an enlarged binding energy in the total system. For 
scalar qnantities, e.g. the effective mass, the separation of phase space has is (in mean held 
approximation) formaly eqnivalent to the introdnction of an additional degree of freedom. 
Using a more realistic interaction derived from the in-medinm Dirac-Brneckner T-matrix this 
effect is present as well, bnt modihed by the snperimposed explicit momentnm dependence 
of the interaction. Dne to the repulsive character of the nuclear forces at high momenta, 
the relative motion of the two cnrrents leads, on the other hand, to a strong additional 
repnlsion in the system. The latter effect originates from the interaction between projectile 
and target nncleons and is taken into acconnt by transport calculation when momentum 
dependent forces are used. The modihcation of the internal mean held of the snbsystems by 
the separation in phase space is, however, not included in standard transport calculations. 
The magnitude of the discussed non-equilibrium ehects is model dependent, however, the 
qnalitative pictnre is model independent. 

In order to discuss the consequences for the eqnation of state in non-eqnilibrium, an 
ehective EOS for colliding matter has been constrncted where the energy of the relative mo¬ 
tion has been snbtracted. It was fonnd as a model independent featnre, that the separation 
of phase space leads to a considerable softening of the EOS in colliding matter. This ehect 
is present over a wide range of densities and streaming velocities. We conclnde that the 
discussed non-equilibrium ehects should be taken into account when transport calculations 
for heavy ion collisions are performed. 
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